*Further Test - Coal Reserve Heterogeneity and Income Level Heterogeneity

clear
eststo clear

*cd "~/Dropbox/Donor Adjustment/Replication/"
use "CombinedData.dta", clear


order WDI_Name id year
xtset id year

*** DATA MANIPULATION

gen log_SO2 = log(SO2+1)
label var log_SO2 "log (SO2)"
replace CO2=0 if CO2 < 0
gen log_CO2 = log(CO2+1)
label var log_CO2 "log (CO2)"
gen log_BOD = log(BOD)
label var log_BOD "log (BOD)"
gen log_gdp = log(gdp)
gen log_PM10 = log(pm10)
label var log_PM10 "log (PM10)"
label var log_gdp "log (GDP)" 
gen log_gdp_pc = log(gdp_pc)
label var log_gdp_pc "log (GDP p.c.)"
gen pop = gdp/gdp_pc
gen log_pop = log(pop)
label var log_pop "log (Population)"
label var polity2 "Polity"
gen time = year-1960
gen time2 = time^2
gen time3 = time^3
gen wto_polity2 = wto * polity2
encode(region), gen(region_id)

gen log_deforest=log(deforest+1)

gen envmin=env_ministry_real
by id (year), sort: replace envmin =. if envmin[_n-1]== 1 
by id (year), sort: replace envmin =. if envmin== 1 & envmin[_n-1]==.

gen FIT1=FIT
by id (year), sort: replace FIT1 =. if FIT1[_n-1]== 1 
by id (year), sort: replace FIT1 =. if FIT1== 1 & FIT1[_n-1]==.
replace FIT1=1 if id==72 & year==1990


gen L1CoalReserve=L1.CoalReserve
label var L1CoalReserve "Coal Reserve (t-1)"

gen middle_income = 0 
replace middle_income = 1 if Classification=="Lower-middle-income economies"
replace middle_income = 1 if Classification=="Lower-middle-income economies"
label var middle_income "Middle-Income Economies"

gen L1middle_income=L1.middle_income
label var L1middle_income "Middle-Income Economies (t-1)"

gen middle_gatt = gatt_wto*middle_income
label var middle_gatt "GATT/WTO*Middle-Income"

gen L1middle_gatt = L1.middle_gatt
label var L1middle_gatt "GATT/WTO*Middle-Income (t-1)"

gen high_income = 0
replace high_income = 1 if Classification=="High-income economies"
label var high_income "High-Income Economies"

gen L1high_income = L1.high_income
label var L1high_income "High-Income Economies (t-1)"

gen high_gatt = gatt_wto*high_income
label var high_gatt "GATT/WTO*High-Income"

gen L1high_gatt = L1.high_gatt
label var L1high_gatt "GATT/WTO*High-Income (t-1)"

gen coal_gatt = gatt_wto*CoalReserve
label var coal_gatt "GATT/WTO*Coal Reserve"

gen L1coal_gatt = L1.coal_gatt
label var L1coal_gatt "GATT/WTO*Coal Reserve (t-1)"

*** DYNAMIC MODEL FOR POLLUTION

gen L1log_SO2=L1.log_SO2
label var L1log_SO2 "log(SO2) (t-1)" 
gen L1log_CO2=L1.log_CO2
label var L1log_CO2 "log(CO2) (t-1)"
gen L1log_BOD=L1.log_BOD
label var L1log_BOD "log(BOD) (t-1)"
gen L1log_PM10=L1.log_PM10
label var L1log_PM10 "log(PM10) (t-1)"

gen L1gatt_wto=L1.gatt_wto
label var L1gatt_wto "GATT/WTO Membership (t-1)"
gen L1log_gdp_pc=L1.log_gdp_pc
label var L1log_gdp_pc "log(GDP p.c) (t-1)"
gen L1log_pop=L1.log_pop
label var L1log_pop "log(population) (t-1)"
gen L1life_expectancy=L1.life_expectancy
label var L1life_expectancy "Life expectancy (t-1)"
label var life_expectancy "Life expectancy"
gen L1polity2=L1.polity2
label var L1polity2 "Polity (t-1)"
gen L1log_environment=L1.log_environment
label var L1log_environment "log(Environmental Aid) (t-1)"
gen L1log_neutral=L1.log_neutral
label var L1log_neutral "log(Neutral Aid) (t-1)"
gen L1log_dirty=L1.log_dirty
label var L1log_dirty "log(Dirty Aid) (t-1)"
label var urban_pop_per "Urban Population $\%$" 
gen L1urban_pop_per=L1.urban_pop_per
label var L1urban_pop_per "Urban Population $\%$ (t-1)"
gen L1electricity_coal_production=L1.electricity_coal_production 
label var L1electricity_coal_production "Electricity from Coal $\%$ (t-1)" 
label var electricity_coal_production "Electricity from Coal"
gen L1electricity_production_renew=L1.electricity_production_renew
label var L1electricity_production_renew "Electricity from Renewable Energy $\%$ (t-1)" 
label var electricity_production_renew "Electricity from Renewables"
gen L1wto=L1.wto
label var L1wto "WTO Membership (t-1)"
label var wto "WTO Membership"

gen coal_wto = wto*CoalReserve
label var coal_wto "WTO*Coal Reserve"

gen L1coal_wto = L1.coal_wto
label var L1coal_wto "WTO*Coal Reserve (t-1)"

gen middle_wto = wto*middle_income
label var middle_wto "WTO*Middle-Income"

gen L1middle_wto = L1.middle_wto
label var L1middle_wto "WTO*Middle-Income (t-1)"

gen high_wto = wto*high_income
label var high_wto "WTO*High-Income"

gen L1high_wto = L1.high_wto
label var L1high_wto "WTO*High-Income (t-1)"


gen L1gatt_wto_informal=L1.gatt_wto_informal
label var L1gatt_wto_informal "GATT/WTO Informal Standing (t-1)"
label var gatt_wto_informal "GATT/WTO Informal Standing"


egen PTA = rowtotal(USPTA JapanPTA EFTAPTA EUPTA)

gen L1auto_emissions_reg=L1.auto_emissions_reg
label var L1auto_emissions_reg "Auto Emission Standards (t-1)"

*cd "~/Dropbox/Donor Adjustment/Replication/check"

****************************
*Coal Reserve Heterogeneity*
****************************

*Table A51: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pol- lution Level
eststo clear
eststo: xi: xtreg log_SO2 gatt_wto coal_gatt   i.year,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto coal_gatt   i.year ,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Pollution1(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{coalpollution1(DID)}) keep(gatt_wto coal_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto coal_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A52: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region, robust 
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region , robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Pollution1(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{coalpollution1(Dynamic)}) keep (L1log_SO2 L1log_CO2 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_SO2 L1log_CO2 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A53: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_BOD gatt_wto coal_gatt   i.year  ,  fe robust
eststo: xi: xtreg log_BOD gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_BOD gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto coal_gatt   i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Pollution2(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{coalpollution2(DID)}) keep(gatt_wto coal_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto coal_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A54: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region , robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region , robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Pollution2(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{coalpollution2(Dynamic)}) keep (L1log_BOD L1log_PM10 L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_BOD L1log_PM10 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps


*Table A55: Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_SO2 wto coal_wto    i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_SO2 wto coal_wto    polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_SO2 wto coal_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto coal_wto      i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto coal_wto    polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto coal_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "PollutionWTO1(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{coalPollutionwto1(DID)}) keep (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A56: Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "PollutionWTO1(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{coalpollutionwto1(Dynamic)}) keep (L1log_SO2 L1log_CO2 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_SO2 L1log_CO2 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A57: Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_BOD wto coal_wto   i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_BOD wto coal_wto   polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_BOD wto coal_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto coal_wto     i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto coal_wto   polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto coal_wto    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "PollutionWTO2(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{coalPollutionwto1(DID)}) keep (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A58: Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "PollutionWTO2(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{coalpollutionwto2(Dynamic)}) keep (L1log_BOD L1log_PM10 L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_BOD L1log_PM10)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A59: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Energy Use
eststo clear
eststo: xi: xtreg electricity_coal_production gatt_wto coal_gatt   i.year  ,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto coal_gatt   i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto coal_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto coal_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Energy(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Energy Use \label{coalenergy(DID)}) keep (gatt_wto coal_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto coal_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A60: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Energy Use
eststo clear
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region , robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1CoalReserve L1coal_gatt  i.year i.region , robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Energy(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Energy Use \label{coalenergy(Dynamic)}) keep (L1electricity_coal_production L1electricity_production_renew L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1CoalReserve L1coal_gatt  L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_coal_production L1electricity_production_renew)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A61: Difference in Difference: Estimating the Heterogenous Effect of WTO on Energy Use after 1995
eststo clear
eststo: xi: xtreg electricity_coal_production wto coal_wto      i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_coal_production wto coal_wto     log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_coal_production wto coal_wto    log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto coal_wto     i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto coal_wto     polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto coal_wto    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "EnergyWTO(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on ENERGY USE after 1995 \label{coalenergywto(DID)}) keep (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto coal_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A62: Dynamic Model: Estimating the Heterogenous Effect of WTO on Energy Use after 1995
eststo clear
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1CoalReserve L1coal_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "EnergyWTO(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Energy Use after 1995 \label{coalenergywto(Dynamic)}) keep (L1electricity_coal_production L1electricity_production_renew L1wto L1CoalReserve L1coal_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_coal_production L1electricity_production_renew)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A63: Probit: Estimating the Heterogenous Effect of GATT/WTO on FIT
set more 1
eststo clear
eststo: xtprobit FIT1 gatt_wto CoalReserve coal_gatt    time time2 time3 i.region_id
eststo:  xtprobit FIT1 gatt_wto CoalReserve coal_gatt    polity2 log_gdp_pc log_pop time time2 time3 i.region_id
eststo:  xtprobit FIT1 gatt_wto CoalReserve coal_gatt    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id
*esttab using "FIT&ENVMIN(probit).tex", booktabs label title(Probit: Estimating the Heterogenous Effect of GATT/WTO on FIT and Environmental Ministry \label{coalFIT&ENV(PROBIT)}) keep (gatt_wto CoalReserve coal_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto CoalReserve coal_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0))  indicate("Time Trend=time" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A64: Probit: Estimating the Heterogenous Effect of WTO on FIT after 1995
eststo clear
eststo:  xtprobit FIT1 wto CoalReserve coal_wto    time time2 time3 i.region_id
eststo:  xtprobit FIT1 wto CoalReserve coal_wto   polity2 log_gdp_pc log_pop time time2 time3 i.region_id
eststo:  xtprobit FIT1 wto CoalReserve coal_wto   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id
*esttab using "FIT&ENVMINWTO(probit).tex", booktabs label title(Probit: Estimating the Heterogenous Effect of WTO on FIT1 & Environmental Ministry after 1995 \label{coalFIT&ENVWTO(PROBIT)}) keep (wto CoalReserve coal_wto  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto CoalReserve coal_wto   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Time Trend=time" )  nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A65: Estimating the Heterogenous Effect of GATT/WTO on Auto Emission Standards
eststo clear
eststo: xi: xtoprobit auto_emissions_reg gatt_wto CoalReserve coal_gatt   polity2 log_gdp_pc L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto CoalReserve coal_gatt   polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto CoalReserve coal_gatt    polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if oecd==0, vce(robust)
*esttab using "AutoEmissionGATT.tex", booktabs label title (Estimating the Heterogenous Effect of GATT/WTO on Auto Emission Standards \label{coalemissiongatt}) keep (gatt_wto CoalReserve coal_gatt  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) order (gatt_wto CoalReserve coal_gatt  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) mtitles ("" "" "Non-OECD" "" "" "Non-OECD") stats (N_g N, labels (Countries Observations) fmt(00)) indicate("Time Trned=time" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A66: Estimating the Heterogenous Effect of WTO on Auto Emission Standards
eststo clear
eststo: xi: xtoprobit auto_emissions_reg wto CoalReserve coal_wto  log_pop polity2 log_gdp_pc time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto CoalReserve coal_wto  log_pop polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto CoalReserve coal_wto  log_pop polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if oecd==0, vce(robust)
*esttab using "AutoEmissionWTO.tex", booktabs label title (Estimating the Heterogenous Effect of WTO on Auto Emission Standards \label{coalemissionwto}) keep (wto CoalReserve coal_wto log_pop polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) order (wto CoalReserve coal_wto log_pop polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) mtitles ("" "" "Non-OECD" "" "" "Non-OECD") stats (N_g N, labels (Countries Observations) fmt(00)) indicate("Time Trned=time" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

****************************
*Income Level Heterogeneity*
****************************

*Table A67: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_SO2 gatt_wto middle_gatt high_gatt   i.year,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto middle_gatt high_gatt   i.year ,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Pollution1(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{econ_pollution1(DID)}) keep(gatt_wto middle_gatt high_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto middle_gatt high_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A68: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region, robust 
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region , robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Pollution1(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{econ_pollution1(Dynamic)}) keep (L1log_SO2 L1log_CO2 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_SO2 L1log_CO2 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A69: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pol- lution Level
eststo clear
eststo: xi: xtreg log_BOD gatt_wto middle_gatt high_gatt   i.year  ,  fe robust
eststo: xi: xtreg log_BOD gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_BOD gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto middle_gatt high_gatt   i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Pollution2(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{econ_pollution2(DID)}) keep(gatt_wto middle_gatt high_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto middle_gatt high_gatt polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A70: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level
eststo clear
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region , robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region , robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Pollution2(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Pollution Level \label{econ_pollution2(Dynamic)}) keep (L1log_BOD L1log_PM10 L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_BOD L1log_PM10 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A71: Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_SO2 wto middle_wto high_wto    i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_SO2 wto middle_wto high_wto    polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_SO2 wto middle_wto high_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto middle_wto high_wto      i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto middle_wto high_wto    polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_CO2 wto middle_wto high_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "PollutionWTO1(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{econ_Pollutionwto1(DID)}) keep (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A72: Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "PollutionWTO1(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{econ_pollutionwto1(Dynamic)}) keep (L1log_SO2 L1log_CO2 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_SO2 L1log_CO2 )  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A73: Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_BOD wto middle_wto high_wto   i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_BOD wto middle_wto high_wto   polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_BOD wto middle_wto high_wto     polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto middle_wto high_wto     i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto middle_wto high_wto   polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg log_PM10 wto middle_wto high_wto    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "PollutionWTO2(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{econ_Pollutionwto1(DID)}) keep (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A74: Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995
eststo clear
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "PollutionWTO2(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Pollution Level after 1995\label{econ_pollutionwto2(Dynamic)}) keep (L1log_BOD L1log_PM10 L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_BOD L1log_PM10)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A75: Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on En- ergy Use
eststo clear
eststo: xi: xtreg electricity_coal_production gatt_wto middle_gatt high_gatt   i.year  ,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto middle_gatt high_gatt   i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop i.year  ,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto middle_gatt high_gatt   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year  ,  fe robust
*esttab using "Energy(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of GATT/WTO on Energy Use \label{econ_energy(DID)}) keep (gatt_wto middle_gatt high_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto middle_gatt high_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A76: Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Energy Use
eststo clear
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region , robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   i.year i.region , robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region, robust
*esttab using "Energy(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of GATT/WTO on Energy Use \label{econ_energy(Dynamic)}) keep (L1electricity_coal_production L1electricity_production_renew L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1middle_income  L1middle_gatt  L1high_income L1high_gatt   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_coal_production L1electricity_production_renew)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A77: Difference in Difference: Estimating the Heterogenous Effect of WTO on ENERGY USE after 1995
set more 1
eststo clear
eststo: xi: xtreg electricity_coal_production wto middle_wto high_wto      i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_coal_production wto middle_wto high_wto    polity2  log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_coal_production wto middle_wto high_wto    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto middle_wto high_wto     i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto middle_wto high_wto     polity2 log_gdp_pc log_pop i.year   if year>=1995,  fe robust
eststo: xi: xtreg electricity_production_renew wto middle_wto high_wto    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per i.year   if year>=1995,  fe robust
*esttab using "EnergyWTO(did).tex", booktabs label title(Difference in Difference: Estimating the Heterogenous Effect of WTO on ENERGY USE after 1995 \label{econ_energywto(DID)}) keep (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto middle_wto high_wto polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A78: Dynamic Model: Estimating the Heterogenous Effect of WTO on Energy Use after 1995
eststo clear
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop  i.year i.region if year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if year>1994, robust
*esttab using "EnergyWTO(dynamic).tex", booktabs label title(Dynamic Model: Estimating the Heterogenous Effect of WTO on Energy Use after 1995 \label{econ_energywto(Dynamic)}) keep (L1electricity_coal_production L1electricity_production_renew L1wto L1middle_income  L1middle_wto  L1high_income L1high_wto   L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1wto L1polity2 L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_coal_production L1electricity_production_renew)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A79: Probit: Estimating the Heterogenous Effect of GATT/WTO on FIT
set more 1
eststo clear
eststo: xtprobit FIT1 gatt_wto middle_income middle_gatt high_income high_gatt    time time2 time3 i.region_id
eststo:  xtprobit FIT1 gatt_wto middle_income middle_gatt high_income high_gatt    polity2 log_gdp_pc log_pop time time2 time3 i.region_id
eststo:  xtprobit FIT1 gatt_wto middle_income middle_gatt high_income high_gatt    polity2 log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id

*esttab using "FIT&ENVMIN(probit)1.tex", booktabs label title(Probit: Estimating the Heterogenous Effect of GATT/WTO on FIT  \label{econ_FIT&ENV(PROBIT)}) keep (gatt_wto middle_income middle_gatt high_income high_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto middle_income middle_gatt high_income high_gatt  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0))  indicate("Time Trend=time" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A80: Probit: Estimating the Heterogenous Effect of WTO on FIT after 1995
eststo clear
eststo:  xtprobit FIT1 wto middle_income middle_wto high_wto high_income    time time2 time3 i.region_id
eststo:  xtprobit FIT1 wto middle_income middle_wto high_wto high_income   polity2 log_gdp_pc log_pop time time2 time3 i.region_id
eststo:  xtprobit FIT1 wto middle_income middle_wto high_wto high_income   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id
*esttab using "FIT&ENVMINWTO(probit)1.tex", booktabs label title(Probit: Estimating the Heterogenous Effect of WTO on FIT after 1995 \label{econ_FIT&ENVWTO(PROBIT)}) keep (wto middle_income middle_wto high_wto high_income  polity2 log_gdp_pc log_pop life_expectancy urban_pop_per) order (wto middle_income middle_wto high_wto high_income   polity2 log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Time Trend=time" )  nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A81: Estimating the Heterogenous Effect of GATT/WTO on Auto Emission Standards
eststo clear
eststo: xi: xtoprobit auto_emissions_reg gatt_wto middle_income middle_gatt high_income high_gatt   polity2 log_gdp_pc L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto middle_income middle_gatt high_income high_gatt   polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto middle_income middle_gatt high_income high_gatt    polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if oecd==0, vce(robust)
*esttab using "AutoEmissionGATT.tex", booktabs label title (Estimating the Heterogenous Effect of GATT/WTO on Auto Emission Standards \label{econ_emissiongatt}) keep (gatt_wto middle_income middle_gatt high_income high_gatt  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) order (gatt_wto middle_income middle_gatt high_income high_gatt  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) mtitles ("" "" "Non-OECD" "" "" "Non-OECD") stats (N_g N, labels (Countries Observations) fmt(00)) indicate("Time Trned=time" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A82: Estimating the Heterogenous Effect of WTO on Auto Emission Standards
eststo clear
eststo: xi: xtoprobit auto_emissions_reg wto middle_income middle_wto high_wto high_income  polity2 log_gdp_pc L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto middle_income middle_wto high_wto high_income  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto middle_income middle_wto high_wto high_income  polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if oecd==0, vce(robust)
*esttab using "AutoEmissionWTO.tex", booktabs label title (Estimating the Heterogenous Effect of WTO on Auto Emission Standards \label{econ_emissionwto}) keep (wto middle_income middle_wto high_wto high_income polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) order (wto middle_income middle_wto high_wto high_income polity2 log_gdp_pc life_expectancy urban_pop_per L1auto_emissions_reg) mtitles ("" "" "Non-OECD" "" "" "Non-OECD") stats (N_g N, labels (Countries Observations) fmt(00)) indicate("Time Trned=time" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps
